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Abstract 

The surface behavior of the pairing gap previously studied for semi-infinite nuclear matter is 
analyzed in the slab geometry. The gap-shape function is calculated in two cases: (a) pairing with 
the Gogny force in a hard-wall potential and (b) pairing with the separable Paris interaction in a 
Saxon- Woods mean-field potential. It is shown that the surface features are preserved in the case 
of slab geometry, being almost independent of the width of the slab. It is also demonstrated that 
the surface enhancement is strengthened as the absolute value of chemical potential \fi\ decreases 
which simulates the approach to the nucleon drip line. 
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1. Introduction 



Recently, the surface behavior of the pairing gap in the ^S'o-channel in semi- infinite 
nuclear matter was investigated independently from two quite different approaches PP, j2]- A 
rather sophisticated approach was used in [T] which starts from the microscopic gap equation 
for semi-infinite nuclear matter with the separable representation [31 E] of the Paris potential 
[S]. The effective pairing interaction Vgg adopted in the gap equation was previously found 
within the Bethe-Goldstone formalism for semi-infinite nuclear matter without any form 
of local approximation [U]. All the calculations were made for two values of the chemical 
potential: fi= — 16 MeV and fi= — 8 MeV. A surface enhancement in the gap A was found, 
the effect being more pronounced for /i = —8 MeV. 

In 12] a more simplified model was used in which nuclear matter was embedded in a semi- 
infinite hard wall potential and the pairing problem was considered in the BCS approximation 
with the Gogny force. Such a simple approach makes it possible to solve the problem to a 
great deal analytically and to examine the coordinate dependence of the pairing gap, pairing 
tensor and correlation energy density in a rather transparent way. A relatively mild surface 
enhancement of all the quantities under consideration was found. As to A, it is of the same 
size as in [T]. 

In this paper we carried out an analogous analysis for a slab of nuclear matter within 
both approaches with the hope that a direct comparison of results can help to clarify the 
general features of the phenomenon under consideration. The slab system is much closer to 
real atomic nuclei than the semi-infinite one and many results can be qualitatively related 
to them. 

The structure of the article is as follows. In Section 2 we extend the model with the 
Gogny force j2] to the case of a hard-wall slab potential. Section 3 contains the extension of 
the model of with the Paris force to the case of slab geometry. The results obtained in 
both models are discussed in Section 4. 

2. Pairing with the Gogny force in the hard wall slab potential. 

Let us consider a slab of nuclear matter embedded in a hard wall potential of thickness 2L 
along the x-direction. Let the origin be at the center of the slab. We start from expanding 
the gap operator in r-space A(r, r') in terms of the wave functions v^k(r) of the hard-wall 
slab potential, 

<^k(r) = y e{L+x) d{L-x) sin K{x-L) e^'^^'"^, (1) 

where k = {/c^jk^}, the quantum number fc^. running over the discrete set of eigenvalues 
kn = 7Tn/{2L), n = 1,2,.... Within the usual BCS approximation the expansion reads 

A(ri, rs) = J2 <^k(ri) <^_k(r2) A(k) , (2) 

k 

where the state | — k) is time-reversed with respect to the state |k). The gap A(k) obeys 
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the BCS equation 



A(p) = - V(p,-p|V^| -k,k) , ^^^^ =, (3) 



where ^(k) = £:(k) — is the single-particle energy relative to the chemical potential and 

(p, -p|l^| - k, k) = JJ dri dv2 v^;(ri) v^lp(r2) ^(ri, ra) </^-k(r2) </^k(ri) (4) 

are the matrix elements of the pairing interaction V{vi,V2)- For the purpose of performing 
calculations analytically we used here, as in [2], the Gogny force Dl 

\/(ri, ra) = j^iW, -B,-H, + e^^'"^-'^)'/"' (5) 

c=l 

with the values of the parameters given in [?]. 

Using the explicit form of the eigenf unctions (P) in Eq. (@)) one obtains 

(p, -p| V^l - k, k) = 4^ (ly, - 5, - + Me) TT e-a^(Pi-fci)/4 

^-^ c=l 

X ^ dxi dx2 9{xi+L) 9{L-xi) 9{x2+L) 9{L-X2) 

X sin(A;a;(xi-L)) sin(A;a;(a;2-Iv)) sin(p^(a;i-L)) sin(p^(a;2--^v)] e~^^^"^^^^''"^ . (6) 

Upon substituting this expression into Eq. Q and, as the s-wave pairing is considered, 
averaging over the angle between vectors p and k, the gap equation can be rewritten as 
follows 



c=l 



" n V 




X 



J J dxidx29{xi+L)9{L-xi)9{x2+L)9{L-X2) 
X sin((A;„(xi-L)) sin(A;„(x2-L)) sin(p„(xi-L)) sin(p^(x2-L)] e^^'^^^''^)'/"' , (7) 



where Io{z) is the modified Bessel function. 

Integrating then over Xi and X2 and introducing the function 



a(v k)- ^v^^^ fc-^^fcV4 



ikar\ „/ikar 2L 



2 a. 



_g2i{p+fc)Lg^£ 



— erf 



ikar 2L 
■ f — 
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ipac 



2L 



- e^^'^P+'^^^erf 



+ 1 erf 



2 a, 

we arrive at the gap equation in the following form 



ipac 2L 



(8) 



1 2 T / 

^ c=l ^ V 

A(A;„, A;_l) 



E 



where 



/e'(A;n,A;^)+A2(A;„,A;j 

Ccb, fc) = ^ Re {s'clp+fc, p+k) + gdp-k, p-k) + gdp+k, -p-k) 
+ 9c{p-k, -p+k) - 2gc{p+k, -p+k) - 2gc{p+k,p-k)} . 



(9) 



(10) 

To reveal spatial behavior of the nonlocal pairing gap operator the Wigner transform of 
the gap is very useful. It reads 



A(R, k) 



dsA(R,s) e 



jks 



'11^ 



where R = (ri+r2)/2 and s = ri— r2. In the case of slab geometry, the Wigner transform 
(jllj) of the gap operator depends only on X which is the component of R perpendicular to 
the surface. Considering A(X, k^, k±) only for X > 0, since it is an even function of X, one 
can easily obtain the following series: 

sin(2A;^(X-L)' 



A(X, k,, k^) = — e{L-X) A(fc„, fcx 



2cos(A;„(X-L) 
sin(2(A;,.-fcJ(X-L)) sin(2(A;,+A;„)(X-L) 



kx 



(12) 



In the bulk, the gap A depends mainly on the total momentum k=\Jk1.+k\. Approx- 
imately, this is true also for the surface region. Within this approximation, we can treat 
the gap at the Fermi surface A{X,kY{X)) as the Wigner transform A(X, /c^., /c_|_) taken at 

k = kpiX), where kp = (37r^p{X)/2) is the local Fermi momentum. We define the 
gap-shape function as the gap at the Fermi surface normalized to one at the center of the 
slab: 

X^(X) = A(x, A;f(X))/A(0, fcF(O)) . (13) 

The function Xp (X) is drawn in Fig. 1 (a and b) for four different values of the half- width 
of the slab: L = 4, 6, 8 and 10 fm. In Fig. la the gap-shape function is calculated for the 
Fermi momentum k-p = 1.4232 fm~^ corresponding to the Fermi energy ep = 42MeV, while 
in calculation of Fig. lb the value of Fermi momentum, kp = 1.4894 fm~^, corresponds to 
Ef = 46 MeV. 
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Figure 1: The gap-shape function Xf(^) 
calculated in the model of pairing with the 
Gogny force in the slab of nuclear matter 
within the hard wall potential for ep = 
42 MeV (panel a) and ep = 46 MeV (panel 
b). The half- width L of the slab is given (in 
fm) by the numbers close to the curves. 



Let us discuss the salient features of the 
shape function Xp(X). The main observation 
is that in the slab, almost independent of its 
thickness, the surface enhancement is not very 
different from the semi-infinite matter case pre- 
viously studied [21 • The enhancement of the 
surface peak is rather moderate, not exceeding 
20% - 30%. In fact it is not clear how much 
of the surface peak comes from Friedel oscilla- 
tions and how much is a genuine enhancement 
of A. From this result and our previous study 
P] we are therefore inclined to conclude that at 
least with the Gogny force pairing is only very 
moderately surface enhanced. As to the depen- 
dence of the gap-shape function on the size of 
the slab, in eq. (fT^ there is a competition be- 
tween two effects. First, with diminishing of the 
size parameter L, the peak value of the numer- 
ator A{X ~ L) becomes higher. Second, the 
denominator (the gap value in the center) also 
gets enhanced. As a result, the L-dependence 
of the gap-shape function turns out to be rather 
smooth. 

Of course, finite nuclei are rather different 
from a slab but qualitatively things should be 
similar. Some results for finite nuclei with the 
Gogny force do exist jHj in the S'n-region but no 
definite conclusions about the surface features 



of A can be drawn from this limited number of 
values. In fact, in finite nuclei like the tin isotopes, we should look at a great number of 
nuclei in the isotopic chain because the behavior of A(r) can fluctuate a great deal passing 
through the open neutron shell in question. 

A detailed study of several long isotopes chains P with density dependent effective 
pairing forces confirmed that, in the framework of a phenomenological approach, it is rather 
difficult to distinguish between two opposite possibilities, the volume pairing and the pairing 
with pronounced surface enhancement. For this purpose, the analysis of some binding energy 
characteristics, such as separation energies, should be accompanied by the study of the 
variations of nuclear radii along the chain. The odd-even staggering phenomenon is especially 
sensitive to the coordinate dependence of A. The microscopic calculation of A should help 
to solve this problem. 

Of course, the conclusion may strongly depend on the employed pairing force and below 
we will investigate the surface behavior of pairing with a separable version of the Paris force. 
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3. Pairing with the Paris force in a slab of nuclear matter 
within a Saxon- Woods potential. 



Now we shall adopt the more realistic one-dimensional Saxon- Woods potential U{X) for 
a slab with the width of 2L symmetrical around the origin X = 0: 



Here Uq is the potential depth in the central region and d is the diffuseness parameter (to 
be more exact, the maximum potential depth is U{0) = Uq/ (1 + 2exp(— L/c/))). The two 
parameters (f/o= — 50 MeV and d=0.65 fm ) are taken to be close to those of real atomic 
nuclei. The half-width parameter L will be changed to examine the size dependence of the 
effect under consideration. 

To avoid a rather cumbersome resolution of the Bogolyubov equations for the nonlocal 
gap [in], as in pp, we use a powerful method [11] (we refer to it as KKC) of solving the 
gap equation for the case of a nonlocal interaction. This method was originally suggested 
for infinite matter where the gap A can be represented as a product A{p) = Apx{p) of 
the constant Ap = A(pf) and the "gap-shape" function x{p) normahzed to x(Pf) = 1- 
Basically the KKC method is a transformation of the gap equation to a set of two coupled 
equations: an integral equation for x(p); which is almost independent of the value of Ap, and 
an algebraic equation for the value Ap. This significantly simplifies the solution of the gap 
equation in infinite matter. In the extension of the KKC method to nonzero temperatures 
fn\ IT^ was used in the case of semi-infinite nuclear matter where the spatial dependence 
of the gap-shape function was also taken into account. Within the KKC method for infinite 
matter, the temperature and momentum parts of the gap function can be factorized as 
follows: A(p, T) = Af(T)x(p). In it was supposed that a similar separation of the 
temperature factor can be made for the semi-infinite system: 



An additional advantage from using the KKC method in this case comes from the possibility 
of finding the normalization factor Ap{T) by solving the gap equation in infinite matter. 

In this paper we use the same ansatz ()15|) for the slab geometry. Unfortunately, in this 
case no direct relation to infinite nuclear matter exists and there is no evident way to find 
the normalization factor without solving the Bogolyubov equations. However, since the gap- 
shape function is of main importance for an analysis of the surface enhancement of the gap, 
we do not calculate here the normalization factor postponing the solution of the Bogolyubov 
equation to a forthcoming publication. 

As far as in the case of the slab geometry all the equations are very similar to those for 
semi-infinite matter, we write down explicitly only those which are necessary for explaining 
our calculations and refer the reader to ^1 E] for details. In symbolic notation, the gap 
equation has the form [121 HI] : 



U{X) 



Uo 



(14) 



1 + exp{{X-L)/d) + exp{-{X+L)/d) ' 



A{xi,X2, k\] T) = Af(T) x(xi, X2; k]_). 



(15) 



A(T) = VP,Ag(r)A(T) 



(16) 
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where Vgg is the effective pairing interaction acting in the model space 5*0 in which the 
superfluid two-particle propagator is defined. 

The separable 3x3 form [3 H] of the Paris potential |3] is used 

V{k,k') = J2K,9r{k')gAk"), (17) 

where k and k' are the relative momenta before and after scattering. The effective interaction 
has a similar separable form which, in notation of ;6], is as follows: 



V^s{xi,X2, X3,X4,kl,k'l;E) = J2 Aij(X, X'; E) gi{k\, x) gj{k 



l2 i\ 

X 



where E=2^, X={xi+X2)/2, x=xi—X2, etc., and gi{k'j_,x) stands for the inverse Fourier 
transform in the x-direction of the form factor gi{k'j^ + k^). The gap-shape factor can be also 
written as 

x{xi,X2;kl) = J2xi{X)gi{kl,x). (19) 

i 

After substituting Eqs. (fT^ . (fTTj) - (fT^ into Eq. (fT^ at T = Tc we obtain the following 
equation for the components Xf 

X.(X) = J2j J dXidX2Au{X,X,;E)BUXuX2,E;T,)xm{X2), (20) 

Im 

where 

dk^ l-Nx,{T)-Nx,{T) 



Bi^{X,X',E;T) = -J2 I 



X 



(27r)2 E - - ex, 

xG'^^^^{kl,X)GZ^^{kl,X'), (21) 

GU(^i'^) = / ynM+^)gi{kl,x)yr.,iX-^)dx, (22) 

NxiT) = (l+exp'-^y\ (23) 

In Eqs. (j2H) - (j2ni) 5 A = (?T,,k^), ex = en+k\/2m, Sn and y„(x) stand for the energies and 
wave functions, respectively, of the 1-dimensional Schrodinger equation with the potential 

m- 

Following the recipe of Ref. [H] for the effective interaction Vgg the propagator em- 
bodies all the single-particle states with negative energies only. Thus, the summation over 
ni,n2 and the integration over k^ are limited by the conditions: e\ < 0, e\i < 0. 

The coefficients Ajj obey the set of integral equations 

Ajj(Xi2, X34; E') = Ajj(5(Xi2— ^34)+ 

+ X! -^^^ j '^^56 Bim {Xi2 , X56 ; E) Amj (X^q , X34 ] E), (24) 



Im 



where Bim are defined by an expression similar to Eq. (PT|) . but without the temperature 
factor, including the states Ai, A2 from the complementary subspace. 

To simplify the calculations we used the local potential approximation (LPA) which has 
turned out to be accurate for semi- infinite nuclear matter [H^ and for nuclear slabs ^H] ■ The 
LPA prescription consists in using for the 2-particle propagator B of the complementary 
space the local momentum approximation, very similar to the Local Density Approximation 
(LDA) [in], for each particle separately: En pl/2m + V{X), En' {p'^)"^ /2m + V{X). This 
type of approach, where the individual particles are treated in semiclassical approximation 
was used in [,17 for examining the response function. This approximation has been shown to 
be very accurate, if one is not interested in fine details. It should be stressed that the LPA is 
only applied to the equation for the effective interaction Vfg, while no local approximation 
is used in the renormalized gap equation (fT^ . That is, the local approximation is used only 
for two-particle states belonging to the complementary space, for which the corresponding 
energy denominators in eq. ()2H) are large. Therefore the individual contribution of each state 
is negligible, and only the sum of a number of such contributions is important. For such a 
sum the semiclassical and local approximations are expected to be accurate. In this respect 
LPA is different from the standard LDA, since in the latter the local approximation is used 
for all two-particle states. 

Within the LPA, the exact values of Bim{Xi,X2;E) are replaced by the set of 
Bl'^it, E;U[X]) calculated for infinite nuclear matter put into the homogeneous potential 
well of the depth Here t=Xi—X2, X={Xi+X2)/2 are the relative and average values 

of the CM coordinates. In the first step of the LPA procedure we calculate the set of vectors 
E = 2/i; Ui) (yi=6V-{i — l)) for a fixed value of the chemical potential /i. In the second 
step, for each value of {X,t), we find Bj^^^Xi, X2) interpolating the values of -B/m(t; Vi) by 
the values of Vi nearest to V{X). After substitution of Bl^^{Xi, X2) into Eq. (j^^ we find 
the LPA prescription Af^^{Xi, X2) for the effective interaction which should be substituted 
into the homogeneous Eq. ()20|) for the gap-shape function. If the critical temperature Tc was 
known, this equation could be solved directly. To find Tc a more general integral equation 
must be considered 



which involves the eigenvalue A(T). Here the abbreviation = Ai?'' is introduced for the 
kernel. The critical temperature can be found from the evident condition A(Tc) = 1. 

The entire calculation scheme is similar to that for semi-infinite matter except for some 
details. First, in the slab case we are dealing with the discrete spectrum e„ in Eq. (j^D). 
Second, due to the obvious reflection symmetry of the slab system in the x-direction, all the 
integral equations under consideration can be readily reduced to a form including positive 
X only. Just as in [Tj, instead of the direct solution of Eq. in the coordinate space, we 
use the Fourier expansion within the symmetrical interval {—Lo,Lq), 




Im 



(25) 



X^{X) = Y: XUniX) 



(26) 



n 



8 



where only the even functions must be retained, /„(a;)= cos(7m(X— Xc)/Lo). The kernels 
Kij{X, X') of Eq. are also expanded in a double Fourier series. Finally we arrive at a 
set of homogeneous linear equations for the coefficients Xi '■ 



3 N 
j=l n'=l 



(27) 



Zp(X): Paris, Woods-Saxon, \i = -8 MeV 



which can be solved by standard numerical methods. Then the components of the gap-shape 
function are found from Eq. ()26j) . 

Instead of analyzing the separate components 
Xi{^) for the separable form (|T7jl of Paris force 
it is more useful to display the local form of the 
gap-shape function which enters the matrix el- 
ements of the gap for states nearby the Fermi 
surface: 

xUX) = j:XiiX)gJk'=kUX)), (28) 




Figure 2: The gap-shape function Xf(^) 
calculated in the model of pairing with the 
Paris force in the slab of nuclear matter 
within the Saxon-Woods potential for fj, = 
-8 MeV (panel a) and n = -4 MeV (panel 
b). The half- width parameter L is given in 
the same way as in Fig. 1. 



where k^iX) = ^2m{fi-U{X)) is the lo- 
cal Fermi momentum {kp{X) = for 
fi—U(X) < 0). We calculated the gap-shape 
function Xp(X) for the same values of the half- 
width of the slab, L = 4, 6, 8 and 10 fm, as in 
the case of the hard wall potential and for two 
values of the chemical potential fi = —8 MeV 
and —4 MeV. As the depth of the Saxon- Woods 
potential was taken Uq = —50 MeV, these two 
values of /i correspond just to the same values of 
the Fermi energy as in the hard- wall case. The 
results are shown in Fig. 2 a and b. 

One observes a surface bump which is much 
more pronounced than in the previous model 
with the Gogny force and the hard wall poten- 
tial. The enhancement is now around 80% - 
100% and, as for the Gogny force, it is quite 
similar to the semi-infinite matter case pp. Two 
values of the chemical potential fi = —8 MeV 
and fi = —4 MeV have been chosen which in ac- 



count of the depth of the Saxon- Woods potential 
of Vo = —50 MeV correspond precisely to the Fermi energies ep = 42 MeV and ep = 46 MeV, 
respectively, of the previous model. Going from fi = —8 MeV to /i = —4 MeV a rather 
important increase of the enhancement of the order of 30% is observed which is much larger 
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than in the case of the hard-wall potential with the Gogny force. On the other hand, there 
is little variation with the thickness of the slab, the maximum of A being a couple of percent 
larger for small slab size. This situation is analogous to the previous model. 

The surface effect in A for the Paris force can be qualitatively explained from the prop- 
erties of the effective interaction Vgg which were analyzed in There it was found that 
Vgff undergoes a sharp variation in the surface region, from almost zero in the bulk to very 
strong attraction in vacuum. In the asymptotic region, the latter coincides with the off-shell 
T-matrix of free A^A^-scattering T{E = 2^) which exhibits a resonant behavior at small 
E. The strong surface attraction and the sharp variation in the surface region are mostly 
responsible for the surface effect of the gap. The //-dependence of the surface effect can be 
explained by the increase of the jump 5Vgg from inside to outside as |//| is decreasing. There 
are two reasons for such an increase. The first one is the fc^-dependence of the form factors 
in Eqs. fll7|) . ()18|) leading to a reduction of V^g in the inner region with increasing values of 
kp. The second one is a pole-like behavior of T[E) at small E which results in an increase 
of V^g with decreasing in the exterior, due to the approach to the virtual pole. One sees 
that both effects work in the same direction resulting in strengthening the surface effect at 
small values of 

4. Discussion, and conclusions. 

In this work we continued our effort to understand the surface behavior of the nuclear 
gap. Previous investigations considered semi-infinite nuclear matter embedded in i) a hard 
wall potential with the Gogny force |2; and ii) a potential of Saxon- Woods shape with the 
separable version of the Paris force [T|. In both cases a surface enhancement was found but 
which is relatively modest in view of what one could expect from LDA. In the present study 
we addressed the question whether finite size effects can strongly alter this situation and 
repeated the former calculation |2] in a slab configuration. 

In the first case of pairing with the Gogny force (Dl) within the slab of nuclear matter, 
the hard wall potential allows to perform most part of calculations analytically. The second 
rather sophisticated case demanded a lot of numerical effort due to the use of the Paris 
interaction and the Saxon- Woods shape of the mean field potential. In both cases a noticeable 
surface effect for the pairing gap was obtained of the same order of magnitude as it was 
previously found in semi-infinite nuclear matter. The shape of the gap in coordinate space 
turned out to be qualitatively similar in both cases, with a significant surface enhancement. 
For the value of the chemical potential fi= — 8 MeV which simulates stable atomic nuclei, 
the enhancement is of the order of 30% for the first model and is almost 100% for the second 
one. A general feature of both models is the rather smooth dependence of the enhancement 
on the slab thickness which is approximately 10% in the first case and only 5% in the 
second one. In both cases, a /i-dependence of the surface effect is found: the enhancement 
coefficient increases as the absolute value of the chemical potential decreases. The latter 
effect is more pronounced for the Paris interaction and Saxon- Woods potential reaching 30% 
with diminishing from 8 MeV to 4 MeV. For the Gogny force and box potential, the 
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corresponding /x-effect is approximately 10%. 

To understand in detail the possible common physical origin of this surface enhancement 
as well as the mentioned differences, it is instructive to consider the gap equation in the 
form (jl6p in which A is expressed in terms of the effective interaction V^g given by Eq. (jl8|) . 
Properties of the effective interaction generated by the Paris force were analyzed in ^H], 
where it was found to undergo a sharp variation in the surface region. It is this sharp 
variation that is mostly responsible for the surface effect in the gap. It is natural that in the 
case of the hard wall box potential the influence of the surface interaction is smaller which 
makes the surface enhancement weaker. 

The /i-dependence of the surface effect is explained by the increase of the jump SV^^ = 
— Vg| with decreasing This increase is caused by two reasons. The first one is the 
strong /c^-dependence of the form factors in Eq. (fT7|) leading to a reduction of for the larger 
values of kp in the bulk. The second one is an increase of V|| with decreasing |/i| caused by 
the pole-like behavior of the T-matrix at small E. Thus, both reasons jointly work towards 
making the surface effect stronger at small values of Qualitatively, the two reasons 
should work also for Gogny force, but the hard wall potential strongly suppresses the second 
reason, diminishing the surface effect itself and its /i-dependence as well. Finally, it should 
be mentioned that a large value of the coherence length of pairing which is comparable with 
the size of the slab makes all the effects under considerations rather smooth. In particular, 
it results in the weak dependence of the surface enhancement of the pairing gap on the slab 
thickness. 

It is worth to point out that the mechanism of the surface enhancement in A considered 
in this paper is different from the one of ref. JHl- The latter is related to a contribution to 
the effective pairing interaction of the virtual exchange by collective surface vibrations. In 
a consistent description of the surface behavior of the pairing gap in nuclei these two effects 
should be considered on equal footing. 

This research was partially supported by Grants No. 00-15-96590 and No. 00-02-17319 
from the Russian Foundation for Basic Research. Two of the authors (E.E.S. and M.V.Z.) 
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